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EXECUTIVE  SUMMARY 

This  report  describes  the  results  of  a  detailed  study  of  the  temporal 
fading  structure  of  scintillating  signals,  particularly  under  conditions 
of  strong  scattering.  In  a  slow- fading  environment,  the  single-bit  error 
rate  depends  only  on  the  probability  density  of  the  signal  fades;  but 
the  fade  coherence  time,  ultimately  determines  the  impact  of  errors 

in  a  data  stream  and  the  effectiveness  (cost)  of  various  coding  schemes 
for  preserving  message  integrity.  In  a  pure  Rayleigh  fading  environment 
the  scintillation  index  is  unity,  and  the  coherence  time  of  the  complex 
signal,  which  is  simply  related  to  T^,  essentially  determines  the  environ¬ 
mental  limitations  to  the  performance  of  any  coherent  or  noncoherent 
signaling  scheme.  Coherence  bandwidth  effects  will  be  discussed  in  a 
separate  report. 

The  principal  result  of  this  study  is  the  verification  of  the  lower 
and  upper  bounds  for  that  are  coming  into  increasing  use  for  system 
performance  evaluation.  Under  conditions  of  weak  scattering,  t  is  a 
monotonic  function  of  the  Fresnel  radius  divided  by  the  effective  scan 
velocity  */z/ve££.  The  details  of  the  functional  dependence  of  Tj  on 
,/z/vef f  depend  on  the  power- law  index  of  the  irregularity  (or  phase) 
spectral  density  function.  In  general,  however,  Tj  <  Vz/veff . 

From  the  Wideband  data  presented  in  this  report  we  find  clear  evidence 
of  the  effects  of  a  varying  spectral  index.  The  wavelength  dependence  of 
Tj  as  determined  from  Wideband  data  recorded  at  Ancon  is  significantly 
different  from  that  determined  from  data  recorded  at  Kwajalein.  This 
difference  is  consistent  with  the  measured  phase  spectral  indices  (Rino 
and  Matthews,  1978)  and  frequency  dependences  of  (Livingston,  1978) 
at  the  two  stations. 

In  spite  of  these  subtleties,  the  data  do  verify  that  <  ./z/v^^. 

The  bound  is  rather  crude,  but  adequate  for  most  engineering  purposes. 
Indeed,  under  conditions  of  weak  scattering,  performance  degradation  is 
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generally  not  severe.  If  a  tighter  bound  is  warranted,  it  can  be 
obtained  by  performing  appropriate  numerical  integrations. 

Under  conditions  of  sufficiently  strong  scattering,  the  fading 
statistics  are  well  approximated  by  the  Rayleigh  model.  In  that  case 
the  structure  of  the  signal  intensity  can  be  characterized  by  a  simple 
asymptotic  limiting  form  of  the  intensity  correlation  function.  Detailed 
analysis  presented  in  Rino  (1978)  shows  that  for  a  phase  spectral  density 
function  of  the  form  „£  »  T^^  This  result  is  strictly 

valid  only  for  1  <  p  <  3.  No  simple  asymptotic  bound  exists  for  p  ^  3. 

To  summarize,  as  one  progresses  from  weak  to  strong  scattering,  t^. 
changes  from  Fresnel  radius  control,  t  <  Jl/v  to  perturbation 
strength  control,  t^.  «  T  ^  .  This  transition  is  clearly  evident  in 

the  Wideband  data. 

However,  a  significant  and  unexpected  result  is  that  under  conditions 
of  strong  scattering  there  is  virtually  a  unique  one-to-one  relation 
between  and  T.  That  is,  there  is  surprisingly  little  dispersion  in  a 
plot  of  Tj.  vs  T.  The  asymptotic  formula  for  the  intensity  correlation 
function  gives  a  tight  lower  bound  for  t^.,  provided  that  p  lies  in  the 
range  of  measured  p  values,  2.3  S  p  £  2.5.  Since  the  calculated  values 
of  decrease  with  increasing  p,  using  p  values  near  3  causes  the 
measured  t^.  value  to  be  severely  underestimated. 

Overall,  the  results  of  our  analysis  further  demonstrate  the 
utility  of  the  phase  screen  model  for  calculating  propagation  effects. 

The  phase  spectral  strength  parameter  T  contains  all  the  necessary 
spatial-to-temporal  conversions  and  the  geometrical  effects  of  propa¬ 
gation  through  highly  anisotropic  media.  It  must  be  kept  in  mind, 
however,  that  the  phase  screen  parameters  are  spatial  averages  over  the 
entire  propagation  path.  The  ramifications  of  this  fact  are  discussed 
in  Section  V. 

The  data  used  in  this  study  were  processed  with  both  a  longer 
detrend  interval  (25  s)  and  a  higher  sample  rate  (250  Hz)  than  are  used 
in  the  routine  summary  processing  of  the  Wideband  data.  This  provided 
an  opportunity  to  double-check  the  results  of  our  routine  data  analysis. 
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The  analysis  does  confirm  that  the  spectral  slopes  are  generally  less 
than  3,  and  that  there  is  a  systematic  difference  between  the  Ancon  and 
Kwajalein  data,  with  the  latter  showing  somewhat  shallower  spectral 
slopes.  Also,  the  longer  detrend  intervals  increased  the  measured 
values  under  strong- scatter  conditions.  This  was  anticipated,  based  on 
the  data  analysis  reported  in  Rino  and  Matthews  (1978). 

For  completeness,  the  phase  scintillation  was  also  analyzed.  The 
simple  linear  dependence  of  the  rms  phase  is  a  good  approximation  to  the 
median  behavior  of  phase  scintillation  for  all  perturbation  levels. 

That  is,  there  is  no  saturation  of  the  phase  scintillation  as  there  is 
with  intensity.  However,  the  dispersion  about  the  mean  increases  steadily 
with  increasing  perturbation  levels.  The  diffraction  effects  in  the 
phase  data  are  evidently  behaving  somewhat  like  additive  noise. 

Finally,  for  engineering  applications,  we  have  computed  the  average 
mean  signal  level  crossing  time,  which  is  easily  measured,  and  compared 
it  with  the  time  delay  to  507c  intensity  decorrelation,  which  is  our 
definition  of  t  .  As  one  should  expect,  the  two  coherence  time  measures 
are  simply  proportional,  particularly  under  strong-scatter  conditions. 
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I  INTRODUCTION  AND  BACKGROUND 


Equatorial  gigahertz  scintillation  is  one  of  a  variety  of  phenomena 
that  are  now  collectively  referred  to  as  "equatorial  spread-F."  The 
discovery  of  gigahertz  scintillation  (Craft  and  Westerlund,  1972)  was 
unexpected,  although  nighttime  equatorial  scintillation  has  been  studied 
for  a  number  of  years  prior  to  Craft  and  Westerlund's  discovery  (Koster, 
1968,  1972). 

The  scintillation  of  transionospheric  radio  signals  can  degrade  the 
performance  of  satellite  communication  links.  Thus,  practical  consider¬ 
ations  have  motivated  a  large  research  effort  over  the  years.  More 
recently,  the  discovery  of  ionospheric  electron  density  depletions 
(Kelley  et  al . ,  1976)  and  their  association  with  equatorial  three-meter 
backscatter  (Woodman  and  La  Hoz,  1976)  have  directed  the  attention  of 
researchers  to  physical  processes  that  can  explain  the  development  of 
the  very  intense  irregularities  that  cause  gigahertz  scintillation 
(Ossakow  and  Chaturvedi,  1978;  Ott,  1978). 

With  regard  to  the  structure  of  the  ionospheric  irregularities  that 
cause  the  scintillation  effects,  it  is  now  known  that  they  encompass  a 
large  continuum  of  scale  sizes  that  can  be  characterized  by  a  power-law 
spectral  density  function  (SDF).  When  the  Fresnel  radius  lies  within 
the  range  of  scale  sizes  characterized  by  the  power-law  continuum,  the 
scintillation  theory  for  both  weak  and  strong  scattering  admits  a  simple 
formulation.  The  weak-scatter  theory  is  described  in  Rino  and  Matthews 
(1978).  The  extension  of  the  theory  to  strong  scattering  is  described 
in  Rino  (1978)  and  the  references  cited  therein. 

In  this  report  we  present  equatorial  data  from  the  Wideband  satel¬ 
lite  that  verify  the  extension  of  the  weak-scatter  theory  to  strong 
scattering.  The  results  are  of  immediate  practical  interest  because 
they  characterize  the  time  structure  of  severely  fading  signals,  which 
is  crucial  in  evaluating  system  performance  and  designing  effective 


mitigants.  Equally  important,  they  provide  a  check  on  the  measured  and 
derived  parameters  that  characterize  the  large  body  of  Wideband  satellite 
equatorial  scintillation  data  that  have  been  accumulated  over  the  past 
two  years. 

Before  presenting  the  data,  however,  we  shall  briefly  review  the 
principal  theoretical  results.  For  a  power- law  irregularity  SDF  of  the 
form 


(K.K  )  =  C  q 
AN  ’  z  sH 
e 


- (2v+l ) 


(1) 


2  2  1/2 

the  scintillation  index  [S^  =  ((i  )  -  (i)  )  /(i)]  under  conditions 

of  weak  scattering  can  be  computed  by  using  the  formula 


=  r^\^(L  sec  @)C  ZV  ^ 
4  e  s 


(**H 


JZ  r/v+0- 


v«  n 


\  2 


(v-0.5) 


it 


(2) 


where 


Xz  sec  4 

4tt 


(3) 


and  3  is  a  geometry- dependent  propagation  factor  [see  Rino  and  Matthews 
1978,  Eq .  (31)].  For  isotropic  irregularities,  ”  =  1;  for  highly 
elongated  irregularities  at  near  normal  incidence,  XT  r(v)/[,./ir  ~(v+F:)]. 
In  Eq.  (2)  rc  is  the  classical  electron  radius,  X  is  the  wavelength,  L 
is  the  layer  thickness,  9  is  tiie  zenith  angle,  and  -  7.  is  the  Fresnel 
radius.  In  Eq.  (3),  zR  is  the  reduced  distance  to  the  equivalent  phase 
screen . 

The  average  structure  of  the  ionospheric  irregularities  is 
characterized  by  the  strength  of  turbulence,  C^,  and  the  spectral  index 
parameter,  v.  In  the  diffraction  theory,  the  admissible  range  of  v 
values  is  0.5  <  ■'  2.5.  If  v  >  2.5,  an  outer-scale  cutoff  must  be 

introduced  to  avoid  divergent  integrals.  Similarly,  if  v  S  0.5,  the 
inner-scale  cutoff  must  be  introduced.  We  note  here  that  Costa  and 
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Kelly  (1976)  have  derived  a  formula  that  retains  an  explicit  dependence  on 
the  outer  scale  when  v  =  1.5.  For  v  =  1.5,  Eq .  (2)  agrees  with  the  Costa- 
Kelly  result  in  the  limit  as  the  outer-scale  wavenumber  approaches  zero. 

In  Figure  1,  which  is  reproduced  from  Rino  and  Matthews  (1978),  we 
show  some  computations  based  on  Eq .  (2)  for  an  overhead  equatorial 
geometry  at  1  GHz.  The  parameter  a  is  the  irregularity  axial  ratio.  For 
a  ^  10  there  is  no  change  in  S^,  and  the  propagation  effects  are  essen¬ 
tially  one-dimensional.  It  can  be  seen  that  for  a  200-km  layer  at  an 

20 

altitude  of  350  km,  Cg  >  10  in  mks  units  is  required  to  produce  signi¬ 
ficant  gigahertz  scintillation. 


FIGURE  1  CALCULATED  GIGAHERTZ  SCINTILLATION  FOR  IDEALIZED 
EQUATORIAL  GEOMETRY 


<> 


The  strength  of  turbulence  is  formally  related  to  the  rms  electron 
2  i 

density,  )  ,  and  the  outer- scale  wavenumber,  q  ,  by  the  normalization 
relation 

c  =  8rt372  <£N2)  q2v"2  r(v4)/r(v-l)  .  (4) 

S  GO 

In  Figure  2,  which  is  also  from  Rino  and  Matthews  (1978),  the  rms  electron 

density  vs  C  for  various  i  =  (2jt/q  )  values  has  been  plotted.  The  rms 
s  oo 

electron  density  values  should  be  interpreted  as  what  would  be  measured 
by  an  in- situ  probe  scanning  along  a  path  of  length  since  in  most 
practical  situations  the  path  length  is  smaller  than  the  outer-scale 
cutoff. 


cs 

FIGURE  2  RMS  ELECTRON  DENSITY  VARIATION  vs.  STRENGTH  OF 
TURBULENCE  FOR  DIFFERENT  VALUES  OF 


/ 
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In  any  case,  it  can  readily  be  seen  that  the  perturbation  levels 
required  to  produce  gigahertz  scintillation  are  very  large  indeed.  It 
should  be  kept  in  mind,  however,  that  insofar  as  the  scintillation  theory 
is  concerned,  Cg  and  v  are  effective  parameters  that  result  from  an 
average  over  the  entire  propagation  path.  Thus,  the  detailed  relationship 
between  Cg,  v,  and  the  corresponding  in- situ  parameters  remains  to  be 
demonstrated . 

Moreover,  to  fully  characterize  the  data,  we  must  take  into  account 
that  only  time-varying  signals  are  actually  measured.  For  a  low-orbiting 
satellite  such  as  Wideband,  the  temporal- to- spatial  conversion  is 
critically  dependent  upon  the  propagation  geometry.  Under  weak-scatter 
conditions,  the  temporal  SDF  of  the  phase  scintillation  has  the  power-law 
form 


cp(f)  =  Tf"P 


(5) 


where 


P  =  2v  (6) 

As  discussed  in  Fremouw  et  al .  (1978),  the  parameters  T  and  p  are  routinely 
measured  in  the  Wideband  data  analysis.  The  relationship  between  Cg  and 
T  [Rino  and  Matthews,  1978],  is 


T 


(L  sec  0)  C  G 
s 


^  r(v) 

(2rt)2v+1  r(v+Sr) 


2v-l. 


eff 


(7) 


where  both  G  and  v^^  depend  on  the  propagation  geometry. 

The  parameter  G  accounts  for  the  scintillation  enhancement  as  the 
propagation  vector  intercepts  the  principal  irregularity  axes,  both 
along  and  (if  appropriate)  across  the  magnetic  field  direction.  The 
parameter,  ve££>  which  has  the  units  of  velocity,  gives  the  temporal- 
to-spatial  coordinate  conversion-- i . e . ,  K  =  (2it£/ve^^) .  For  isotropic 
irregularities  or  cross-field  scan  directions,  =  vj_,  where  v±  is 

the  component  of  the  scan  velocity  perpendicular  to  the  propagation 
direction.  In  general,  v^^  ^  vx* 
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As  Cg  increases,  the  measured  value  of  ultimately  falls  below 
that  predicted  by  Eq.  (2).  Numerical  calculations  for  a  strict  power-law 
environment  predict  values  increasing  to  a  maximum  value  somewhat 
larger  than  unity  and  then  converging  to  unity  from  above  [Rumsey  (1975); 
Marians  (1975)3-  For  Rayleigh  fading,  =  1,  although  =  1  is  not  a 
sufficient  condition  to  guarantee  Rayleigh  fading.  When  >  1,  there  is 
an  excess  of  signal  enhancements  over  that  expected  for  pure  Rayleigh 
fading.  The  degree  to  which  this  effect  occurs,  depends  on  the  spectral 
index  v.  Spectra  with  large  v  values  intensify  the  scale-size  regime 
near  the  Fresnel  frequency,  which  produces  so  called  strong  focusing. 

When  this  occurs,  the  index  exceeds  unity. 

Pure  Rayleigh  fading,  however,  is  the  most  commonly  used  engineering 

model.  If  x  and  y  are  zero-mean  and  uncorrelated  but  identical  Gaussian 

processes,  then  the  complex  signal  v  =  x  +  iy  has  a  Rayleigh  distribution 

2 

for  its  intensity,  I,  (I  =  |vj  )  and  a  uniformly  distributed  phase  between 
0  and  2rt  radians.  The  quantity  (w  )  is  called  the  mutual-coherence 
function,  and  the  relation 

<Il'>  -  <I>2  =  |<vv'*x|2  =  4<xx')2  (8) 

2 

is  easily  verified.  Since  <l)  =  2<x  >  =  1,  it  follows  from  the  definition 
of  and  Eq.  (8)  that  =  1  as  noted  above. 

It  is  also  well  known  (Fejer,  1953;  Bramley,  1954)  that  for  a 
multiply  scattered  wavefield  in  a  uniform  medium, 

(w'">  =  exp  j-  i-D(Ap)}  (9) 

where  D(^p)  is  the  phase  structure  Sanction 

D(^)  =  <(6<2  -  60 '  )2v  ,  (10) 

and 

50  =>  rg\  (  •ANe  dl  (11) 
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The  integral  is  evaluated  along  the  propagation  path.  Thus,  under 
conditions  of  sufficiently  strong  scattering  where  ~  1  and  a  relation¬ 
ship  such  as  Eq.  (8)  is  expected  to  hold,  Eq.  (9)  can  be  used  to  derive 
a  simple  asymptotic  expression  for  the  fade  coherence  time. 

In  Rino  (1978)  a  detailed  analysis  of  the  intensity  correlation 
function  under  conditions  of  strong  scattering  is  presented.  The  analysis 
shows  that  in  a  strict  power-law  environment,  the  asymptotic  strong- 

—4 

scatter  limit,  Eq.  (8),  is  valid  as  long  as  V  <  1.5.  If  v  >  1.5,  D(/r  ) 
does  not  exist  in  the  limit  of  an  arbitrarily  large  outer  scale,  and  the 
intensity  correlation  function  as  well  as  the  mutual-coherence  function 
retain  dependences  on  the  Fresnel  radius  no  matter  how  strong  the  pertur¬ 
bation  becomes.  This  is  not  a  contradiction  to  the  Fe jer-Bramley  result 
because  a  scattering  medium  with  an  arbitrarily  large  outer-scale  para¬ 
meter  is  not  strictly  uniform. 

Nonetheless,  in  the  simple  case  where  v  1.5,  the  time  structure 
of  the  fading  under  strong-scatter  conditions  asymptotically  approaches 


<Il'>  =  exp  |-D(vs5t)}  . 


(12) 


It  is  shown  in  Rino  (1978)  that 


D(vg5t)  =  T 


4rt2v~ °* 5  r(1.5-v) 


(2v-d  rco 


5t 


2V-1 


(13) 


Thus,  the  intensity  coherence  time  under  conditions  of  strong  scattering 
admits  a  simple  characterization  that  can  be  used  to  verify  the  measured 
T  and  v  =  p/2  values.  We  shall  see  that  the  results  are  in  good  agree¬ 
ment  with  the  Wideband  satellite  data. 


II  THE  DATA  BASE 


To  obtain  a  representative  data  base  of  intense  scintillation,  we 
have  selected  sets  of  the  most  severely  disturbed  passes  recorded  at  the 
Ancon  (11°46'S,  77°09'W)  and  Kwajalein  (9°24'N,  167°28'E)  tracking 
stations.  The  specific  passes  are  listed  in  the  Appendix.  The  occurrence 
statistics  of  equatorial  scintillation  are  discussed  in  Livingston  (1978). 

In  the  routine  summary  processing  that  is  applied  to  all  the  Wide¬ 
band  satellite  data,  a  10-s  detrend  filter  is  used  with  a  100-Hz  sample 
rate  (Fremouw  et  al.,  1978).  To  fully  accoimodate  the  spectral  broadening 
associated  with  intense  scintillation,  however,  it  is  desirable  to 
increase  the  temporal  resolution.  Moreover,  as  discussed  in  Section  IV 
of  Rino  and  Matthews  (1978),  the  10-s  detrend  interval  is  evidently  too 
short  to  fully  accommodate  the  low-frequency  content  of  the  most  severe 
equatorial  intensity  scintillation. 

For  this  study,  therefore,  the  L-band,  center  UHF,  and  VHF  channels 
have  been  reprocessed  at  a  250-Hz  sample  rate  using  a  25-s  detrend 
filter.  The  phase  data  were  then  spectrally  analyzed  to  extract  the  T 
and  p  parameters  [see  Eq .  (5)]  as  described  in  Fremouw  et  al .  (1978). 

To  characterize  the  time  structure  of  the  signal,  the  time  lag  to  507 
intensity  decorrelation  was  measured  as  well  as  the  mean  intensity 
crossing  rate.  The  summary  parameters  were  measured  on  non-overlapping 
16-s  time  intervals. 

To  illustrate  the  effects  of  the  detrender  cutoff,  in  Figure  3  a 
segment  of  UHF  data  processed  using  both  the  standard  10-s  detrend  filter 
and  a  25-s  detrend  filter  is  shown.  To  the  eye,  the  intensity  data  are 
identical.  Thus,  there  is  very  little  spectral  content  in  the  intensity 
data  below  0.1  Hz.  This  intrinsic  high-pass  structure  of  the  intensity 
scintillation  is  caused  by  the  diffraction  process,  which  strongly 
suppresses  Fourier  modes  with  spatial  periods  longer  than  the  Fresnel 
radius . 
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FIGURE  3  DATA  SEGMENT  SHOWING  EFFECTS  OF  CHANGING  DETREND  FILTER  CUTOFF 


Since  the  phase  data  essentially  map  the  integrated  electron  content 
along  the  propagation  path,  however,  the  low-frequency  modes,  which  are 
the  most  intense  in  a  power-law  environment,  dominate  the  phase  structure. 
Thus,  there  is  a  dramatic  increase  in  the  low-frequency  content  of  the 
phase  scintillation  data  when  the  detrend  interval  is  increased. 

This  intrinsic  nonstationarity  in  the  phase  data  is  easily 
characterized  in  terms  of  the  T  and  p  parameters.  For  example,  the 
measured  phase  variance  is 

CO 

<5d>2)  =  2T  f  f'P  df 


=  ^TTcP’1  (14) 

P-1  c 

where  tc  is  the  detrend  interval.  The  Wiener  process,  which  is  used  to 
model  Brownian  motion,  has  a  similar  property  [Uhlenbeck  and  Ornstein, 
19301, 

Unfortunately,  there  is  no  theory  that  either  adequately  accommodates 
diffraction  effects  in  phase  scintillation  or  even  addressed  the  phase 
structure  under  conditions  of  strong  scattering.  In  our  analysis,  only 
the  zeroth-order  approximation,  under  which  diffraction  effects  are 
entirely  neglected,  has  been  considered.  Thus,  the  temporal  phase  SDF 
has  the  form  Tf  P;  p  is  independent  of  wavelength  and  T  varies  quadrati- 
cally  with  wavelength. 

It  follows  from  Eq .  (14)  that 

2  1/2 

-  =  (^  )  *  A.  .  (15) 


To  correct  for  a  finite  reference  frequency,  Eq.  (15)  must  be  multiplied 
2 

by  1  -  (f/f^)  ,  where  f  =  c/X,  and  f  ^  is  the  reference  frequency.  Thus, 
the  wavelength  dependence  of  the  rms  phase  gives  one  simple  test  of  the 
adequacy  of  the  simplest  theory.  In  Figure  4  we  show  scatter  diagrams  of 


rms  phase  for  L-band  vs  UHF  and  UHF  vs  VHF  from  Che  Ancon  data.  The 
corresponding  scatter  diagrams  for  the  Kwajalein  data  are  shown  in 
Figure  5. 

When  the  rms  phase  is  less  than  a  few  radians,  the  data  adhere  very 
closely  to  the  theoretical  curves.  For  larger  perturbation  levels  there 
is  steadily  increasing  but  essentially  uniform  scatter  about  the 
theoretically  predicted  means.  There  is  evidently  no  saturation  effect 
as  there  is  with  the  intensity  scintillation  index. 

Intuitively,  it  is  clear  that  diffraction- induced  departures  from 
the  simple  theory  are  caused  by  the  rapid  phase  changes  that  accompany 
deep  fades.  Thus,  it  is  mainly  the  high  frequencies  that  are  affected, 
and  the  increasing  scatter  in  Figures  4  and  5  for  the  larger  perturbations 
is  essentially  additive  "diffraction"  noise. 

It  should  be  kept  in  mind  that  there  is  no  significance  to  the 

2  2 

absolute  rms  phase.  From  Eq .  (14),  however,  for  p  =  3,  )  =  Ttc- 

Thus,  the  one  radian  level  where  diffraction  effects  become  significant 

_  3 

corresponds  to  T  =  1.6  X  10  ,  or  -28  dB.  It  is  convenient  to  express  T 

in  dB  because  it  is  numerically  equal  to  the  phase  SDF  evaluated  at  1  Hz. 

In  Figures  6  and  7  the  averages  of  the  measured  p  indices  are  shown 
for  the  Ancon  and  Kwajalein  data  sets,  respectively.  The  data  have  been 
ordered  by  scintillation  level  to  show  any  systematic  changes  with 
perturbation  strength.  For  very  weak  scintillation  (S^  -  0.1)  the 
system  noise  obscures  the  phase  SDF  and  p  is  artificially  decreased. 

For  large  perturbation  levels  (S^  **  0.8),  diffraction  effects  act  to 
reduce  the  phase- spec tral  slope. 

The  shaded  regions  in  Figure  6  denote  the  standard  deviation  of  the 
UHF  data.  The  data  in  Figures  6  and  7  show  that  for  comparable  pertur¬ 
bation  levels,  the  measured  p  indices  are  indeed  wave! ength- invariant . 
Moreover,  the  median  p  values  for  the  intermediate  range  where  neither 
noise  contamination  nor  diffraction  significantly  affect  the  measurements, 
are  consistently  lower  than  p  =  3.  Finally,  we  note  that  there  is  a 
small  but  systematic  difference  between  the  spectral  slopes  from  the 
Ancon  data  and  the  Kwajalein  data.  For  the  former,  p  ~  2,6,  and  p  ~  2.4 
for  the  latter. 
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FIGURE  5  SCATTER  PLOTS  OF  o ^  FOR  KWAJALEIN  DATA 


FIGURE  6  AVERAGE  OF  MEASURED  ANCON  p  INDICES  AT  VHF,  UHF,  AND  L-BAND 
PLOTTED  AGAINST  S4 

The  results  confirm  the  earlier  findings  reported  in  Rino  and 
Matthews  (1978)  that  were  based  on  the  standard  summary  data.  Moreover, 
Livingston  (1978)  reports  a  difference  in  the  wavelength  dependence  of 
the  average  intensity  scintillation  measured  at  Ancon  and  Kwajalein. 

For  roughly  comparable  L-band  scintillation  levels,  there  is  significantly 
more  VHF  scintillation  observed  at  Ancon  than  at  Kwajalein.  This 
difference  can  be  fully  accounted  for  by  the  slight  difference  in  spectral 
slopes,  with  the  steeper-slope  spectra  producing  more  intensity  scintil¬ 
lation  for  the  same  perturbation  level  (see  Figure  1). 

Turning  now  to  the  intensity  scintillation,  in  Figures  8  and  9  the 
two- frequency  scatter  diagrams  of  measured  values  that  correspond  to 
Figures  4  and  5,  respectively,  are  shown.  Under  conditions  of  weak 
scatter  (say  0.4)  there  is  a  sensibly  linear  relation  between  the 

values  measured  at  two  different  frequencies.  This  is  in  agreement 
with  the  weak-scatter  theory,  which  predicts 
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FIGURE  7  AVERAGE  OF  MEASURED  KWAJALEIN  p  INDICES  AT  VHF,  UHF,  AND  L-BAND  PLOTTED 
AGAINST  S4 


£±2 

S  CC  \  4 
4 


(16) 


[see  Eq.  (2)  and  recall  that  p  =  2v ] . 

Unfortunately,  the  scatter  in  the  measured  data  points  is  such  that 
the  departure  of  the  spectral  slopes  from  p  =  3  and  the  subtle  differences 
between  Ancon  and  Kwajalein  are  not  immediately  apparent.  On  the  other 
hand,  if  Figures  8(a)  and  9(a)  are  compared,  it  can  be  seen  that  there 
are  relatively  fewer  Kwajalein  data  points  in  the  intermediate  range. 

As  discussed  in  Livingston  (1978),  this  effect  is  consistent  with  the 
small  difference  in  spectral  slopes. 

Under  strong- sc  at  ter  conditions  there  is  essentially  uniform  scatter 
about  S,  =  1 .  We  believe  that  these  are  statistical  fluctuations  rather 
than  evidence  of  the  strong  focusing  effect  discussed  in  Section  I. 

Indeed,  in  a  one-dimensional  medium  the  strong  focusing  effect  does  not 
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FIGURE  9  SCATTER  PLOTS  OF  S a  FOR  KWAJALEIN  DATA 


occur  unless  p  &  3.  In  any  case,  these  subtleties  were  suppressed  in 
the  standard  summary  data  reported  in  Rino  and  Matthews  (1978)  because 
of  the  10-s  detrend  filter  cutoff. 

To  summarize,  a  representative  set  of  disturbed  equatorial  passes 
has  been  processed  to  obtain  both  high  time  and  frequency  resolution. 

For  small  perturbation  levels,  the  first-order  phase  and  amplitude  signal 
moments  follow  the  weak-scatter  theory.  As  the  perturbations  increase, 
the  rms  phase  continues  to  follow  the  first-order  theory,  albeit  with 
considerably  more  scatter.  The  values,  on  the  other  hand,  saturate 
about  a  median  value  of  =1.  There  is  no  evidence  of  strong  focusing, 
which  is  in  agreement  with  the  theory  and  further  confirmation  of  the 
shallow  p  indices. 

The  structure  of  the  equivalent  phase-changing  screen  is  completely 
characterized  by  the  phase  SDF  through  the  T  and  p  parameters.  The  T 
parameter  combines  the  effects  of  absolute  turbulent  strength,  and  the 
conversion  of  spatial  to  temporal  structure  in  an  anisotropic  medium. 

The  diffraction  effects  are  characterized  by  the  Fresno  parameter,  Z. 
Because  of  the  very  large  extent  of  the  power- lav  continuum,  no  signi¬ 
ficant  error  is  made  in  neglecting  both  the  inner-  and  outer-scale 
spectral  cutoffs.  That  is,  these  cutoffs  are  well  outside  the  range  of 
scale  sizes  that  exert  a  measurable  effect  on  the  scintillation. 
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Ill  THE  TEMPORAL  STRUCTURE  OF  INTENSITY  SCINTILLATION 


To  characterize  the  temporal  structure  of  intensity  scintillation, 
the  time  delay  to  50 °L  decorrelation,  t^,  has  been  used.  Formally, 

<1(tw  + 1 , 0,5  <17, 

s4 

where  it  is  understood  that  (i)  =  1.  For  some  applications,  however,  it 
is  desirable  to  have  a  simpler  measure  of  the  fade  coherence  time.  Thus,' 
the  reciprocal  of  the  upward  (or  downward)  mean  level  crossing  rate,  t 
has  also  been  computed. 

One  expects  a  simple  monotonic  relation  between  t^.  and  ,  with 
generally  greater  than  t  .  Consider  that  a  series  of  positive  and 
negative  impulses  has  a  very  short  coherence  time,  but  necessarily  a 
longer  time  between  impulses.  To  verify  this  behavior,  in  Figure  10 
scatter  plots  of  vs  r  are  shown  for  the  Ancon  data.  The  corresponding 
Kwajalein  data  are  shown  in  Figure  11.  Most  of  the  points  follow  an 
essentially  linear  relation.  The  dispersion  is  particularly  small  for 
the  stronger-scatter  cases  where  the  coherence  times  are  small. 

There  are,  however,  a  large  number  of  data  points  with  v>  Tc* 

These  points  are  confined  to  conditions  of  weaker  scatter  where  the 
fluctuations  about  the  mean  intensity  level  are  comparatively  small.  To 
verify  this,  the  Kwajalein  data  shown  in  Figure  10  have  been  replotted 
(Figure  12)  with  all  data  points  for  which  S  0.8  removed.  The  few 
remaining  anomalies  are  attributed  to  nonstationar ities  which  invariably 
occur  as  the  propagation  path  passes  through  regions  of  varying  pertur¬ 
bation  level . 

Now,  under  conditions  of  weak  scatter,  the  intensify  correlation 
function  is  given  by  the  formula 
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FIGURE  11  SCATTER  PLOT  OF  tc  vs.  Tj  FOR  KWAJALEIN  DATA 


in  terms  of  the  actual  scan  velocity,  the  propagation  geometry,  and  the 
irregularity  anisotropy  (Rino  and  Livingston,  1978,  Section  II). 

Unfortunately,  rather  messy  numerical  integrations  would  be  required 

to  evaluate  Eq.  (18).  Nonetheless,  it  can  be  seen  that  t  is  independent 

of  perturbation  strength.  Moreover,  one  expects  to  be  a  monotomic 

function  of  Jt. / v  where  v  is  the  effective  scan  velocity  discussed 

eff  eff 

in  Section  1.  To  verify  this  dependence,  in  Figure  13  a  scatter  diagram 
of  login  t  vs  ^/v  .  is  shown. 
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FIGURE  12  REPLOT  OF  DATA  FROM  FIGURE  10  WITH  DATA  POINTS  FOR 
S4  <  0.8  REMOVED 


The  Fresnel  radius  ■  Z  was  calculated  for  a  350-km  reference  altitude. 
The  effective  scan  velocity  was  calculated  for  50:1  rod-like  irregularities. 
Superimposed  on  the  scatter  diagram  is  the  theoretical  curve  =  K.  Z/v^^ 
for  K  =  0.16.  The  value  of  K,  however,  is  model-dependent.  Thus,  the 
best  one  can  do  for  predictive  purposes  is  to  take  ■.  Z/v  ^  as  an  upper 
bound  to  the  fade  coherence  time  under  conditions  of  weak  scatter. 

Insofar  as  the  t  vs  . Vz/v  behavior  is  concerned,  the  scatter 

I  eft 

diagrams  for  Kwajalein  and  Ancon  data  show  no  perceptible  difference. 

Thus,  we  have  shown  only  the  Kwajalein  data.  On  the  other  hand,  if  the 
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FIGURE  13  SCATTER  DIAGRAM  OF  LOGiq  T|  vs.  J~ Z/ve^.  Superimposed  is 
a  plot  of  7|  =  K  V7/veff. 


Tj  values  measured  at  each  frequency  are  plotted  separately,  differences 
between  Ancon  and  Kwajalein  data  that  are  consistent  with  the  different 
spectral  slopes  do  show  up. 

In  Figure  14,  t  vs  plots  are  superimposed  for  the  Ancon  VHF, 

UHF,  and  L-band  data.  For  weak  scatter--say,  i  0.4--t^.  is  independent 
of  but  varies  monotonical ly  with  wavelength  because  of  the  dependence 
on  the  Fresnel  radius.  As  the  perturbation  strength  increases,  the 
wavelength  dependence  disappears,  giving  way  to  a  simple  dependence  on 
perturbation  strength,  as  will  be  shown  in  more  detail  shortly.  This 
general  behavior,  however,  is  clearly  seen  in  Figure  1A. 


AVERAGE  T\  VALUES  FOR  ANCON  DATA  vs. 


The  corresponding  data  for  Kwajalein  are  shown  in  Figure  15.  The 
general  behavior  is  similar  to  that  of  the  Ancon  data.  However,  the 
weak- scatter  wavelength  dependence  is  much  less  pronounced.  Indeed, 
it  is  virtually  lost  in  the  fluctuations  of  the  measurements  themselves. 
This  behavior  is  attributable  to  the  generally  flatter  phase  SDFs 
observed  at  Kwajalein,  as  discussed  in  Section  II. 

Let  us  now  consider  the  intensity  coherence  time  under  conditions 
of  strong  scattering.  Taking  Eqs.  (12)  and  (13)  as  guidelines,  it 
follows  that 


where 


(20) 


for  507„  decorrelation.  Thus,  log^  should  vary  linearly  with  log.^  T, 
with  a  slope  that  depends  on  the  spectral  slope  Y  =  p/2. 

Thus,  in  Figures  16  and  17,  log-log  plots  of  t^.  vs  T  are  shown  for 
the  Ancon  and  Kwajalein  data,  respectively.  Only  points  for  >  0.8 
have  been  included.  The  first  feature  to  note  in  both  plots  is  the  very 
small  amount  of  dispersion  in  the  data  for  the  larger  perturbation 
levels.  There  is  effectively  a  one-to-one  relation  between  T  and 

under  strong-scatter  conditions.  Moreover,  the  relationship  follows 

£ 

the  power- law  form  t  «  T  with  an  index  e  =  2v  -  1  very  close  to  the 
value  inferred  from  the  measured  p  indices. 

In  Figures  16  and  17,  theoretical  curves  based  on  Eq.  (19)  are  also 
shown  for  p  =  2.9  and  p  equal  to  the  index  corresponding  to  the  least- 
squares  fit  to  the  respective  data  sets.  It  can  be  seen  that  the  p  =  2.9 
curves  perdict  much  shorter  coherence  times  for  a  given  T  level  then 
are  actually  observed.  The  curves  corresponding  more  closely  to  the 
measured  p  indices  still  predict  somewhat  smaller  coherence  times  than 
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IGURE  15  SEMI-LOG  PLOT  OF  AVERAGE  r,  VALUES  FOR  KWAJALEIN  DATA  vs. 


FIGURE  16  LOG-LOG  PLOT  OF  r,  vs  T  FOR  ANCON  DATA  TOGETHF  R  WITH  LEAST-SOU ARFS 
FIT  TO  DATA  AND  THEORETICAL  CURVES  DERIVtD  FROM  Eq  1191 


are  actually  observed.  Moreover,  the  (inferences  between  Ancon  and 
Kwajalein  as  i ar  as  the  inferred  p  index  values  are  verv  snail. 

These  small  variances  between  the  theory  and  the  data  are  attri¬ 
butable  to  the  I  act  that  l>|.  (191  is  an  asymptotic  result  that  applies 

strictly  on  1 v  in  the  limit  .is  the  perturbation  becomes  arbitrarily  large. 
Moreover,  the  convergence  depends  on  the  spectral  index,  with  the  smaller 
indices  converging  more  1  •  1  y .  Overall,  tin  internal  consistency  of 

the  data  and  its  agreement  with  the  theory  are  very  good  indeed. 

To  summarize  the  practical  ramifications  ot  this  section,  we  have 


demonstrated  two  bounds  for 


I'nder  conditions  of  weak  scattering 
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FIGURE  17  LOG-LOG  PLOT  OF  T |  vs.  T  FOR  KWAJALEIN  DATA  TOGETHER  WITH  LEAST-SQUARES 
FIT  TO  DATA  AND  THEORETICAL  CURVES  DERIVED  FROM  Eq  (19) 


rj  is  independent  of  perturbation  strength  and  bounded  from  above  by 
Z/v  This  bound,  however,  is  rather  crude.  Under  conditions  of 

strong  scattering  a  much  tighter  lower  bound  on  tj  is  given  by  Eq .  (19). 
Thus,  for  predictive  purposes  it  is  reasonable  to  take 


Fortunately,  a  tight  bound  on  tj  under  weak-scatter  conditions  is 
largely  unnecessary  because  the  system  effects  are  not  severe. 


IV  DISCUSSION 


The  data  analysis  presented  in  this  report  demonstrates  the  utility 
of  the  phase- screen  model  for  interpreting  ionospheric  scintillation. 

The  results  are  particularly  significant  in  light  of  new  theoretical 
results  that  extend  the  theoretical  results  based  on  the  phase-screen 
model  into  the  strong-scatter  regime.  This  regime,  which  was  previously 
thought  to  be  mathematically  intractable,  has  yielded  results  that 
reproduce  the  data  much  more  faithfully  than  the  corresponding  weak- 
scatter  results. 

Thus,  in  the  weak-scatter  regime  where  the  intensity  coherence  time 
is  controlled  by  fz/v eff,  the  data  show  both  a  high  degree  of  dispersion 
and  a  sensitivity  to  the  precise  value  of  the  spectral  index  that  is  very 
difficult  to  model.  By  comparison,  under  strong-scatter  conditions  only 
the  phase  turbulent  strength,  T,  and  the  spectral  index,  p,  need  be  known 
to  completely  specify  the  intensity  autocorrelation  function.  Moreover, 
the  formula  is  a  simple  one,  and  it  accurately  reproduces  the  measured 
coherence  times. 

The  broader  question  of  how  the  phase-screen  parameters  relate  to 
the  actual  in-situ  irregularity  structures  deserves  some  careful  consi¬ 
deration.  The  interpretation  of  scintillation  has  been  strongly  in¬ 
fluenced  by  the  work  of  Dyson  et  al.  (1974)  who  reported  one-dimensional 
in-situ  SDFs  witli  power- law  indices  somewhat  less  than,  but  very  near, 

2  (1.84  ±-'-0.1).  Most  researchers  have  been  content  to  use  the  nearest 
integral  value  2  for  both  theoretical  calculations  and  data  interpretation. 

Moreover,  it  is  well  demonstrated  that  both  the  gradient-drift  and 
the  Rayleigh-Taylor  (gravitational)  instability  generate  steep  gradients. 

If  such  gradients  are  idealized  by  actual  discontinuities,  then  the  high- 
frequency  portion  of  the  one-dimensional  SDF  will  certainly  have  the 
form  (Costa  and  Kelley,  1976). 
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To  illustrate  this,  consider  a  collection  of  randomly  distributed 
rods  on  a  line.  If  the  diameters  of  the  rods,  i,  are  distributed 
according  to  the  probability  law  p(£),  then  it  is  easily  shown  that  the 
one-dimensional  SDF  has  the  form 


cc 

9 

• 

pa> 

sin  K  1 

X 

=  K-2 

X 

1-2  cos  (2K^)  pU>  di 

CM 

Jo 

L  x  j 

•4) 

(22) 


For  any  smooth  function  p(i),  the  second  term  will  decay  more  rapidly 
than  Ky2,  Hence,  for  sufficiently  large  K^,  ^(K^)  }Q2. 

On  the  other  hand,  one  can  easily  generate  a  reasonable  size  distri¬ 
bution  that  will  produce  a  broad  spectral  regime  with  a  spectral  index 
between  0  and  2.  A  more  realistic  distribution  of  randomly  distributed 
rods  on  a  plane  has  been  treated  by  Chesnut  (1978)  with  similar  results. 

If  such  a  simplified  model  is  indeed  valid,  the  important  phenom- 

_2 

enology  question  is,  what  spatial  wavenumbers  encompass  the  and  the 

“  2-f-s 

1^  regimes?  The  simplest  relationship  between  the  in- situ  and  phase 
SDFs  (Rino  and  Matthews,  1978)  gives  a  p  index  one  larger  than  the 

spectral  index  corresponding  to  the  one-dimensional  in-situ  SDF.  Thus, 

-2 

Kj.  corresponds  to  p  =  3. 

Since  the  Wideband  satellite  data  have  consistently  shown  p  values 

less  than  3  for  both  the  auroral  zone  and  equatorial  data,  it  is 

tempting  to  conclude  that  the  spatial  wavenumbers  that  cause  scintil- 
-  2+e 

lation  lie  in  the  regime.  Stated  another  way,  the  Wideband  data 

_  o 

suggest  that  the  portion  of  the  in-situ  SDF  cannot  persist  below 

spatial  wavenumbers  corresponding  to  ~  2 00  m. 

However,  Wittwer  (1978)  has  pointed  out  that  both  actual  source 
variations  and  geometrical  smearing  duo  to  the  varying  magnetic  aspect 
angle  along  the  propagation  path  can  upset  the  simple  relation  between 
p  and  the  one-dimensional  in-situ  index.  Studies  are  currently  underway 
to  assess  these  effects  both  bv  looking  at  in-situ  data  and  performing 
theoretical  computations.  Thus,  a  definitive  link  must  await  the 
culmination  of  these  efforts. 


5  b 


Nonetheless,  the  fact  that  verifiable  systematic  differences  in 
structure  are  present  in  the  Ancon  and  Kwajalein  data,  which  have 
virtually  identical  propagation  geometries,  strongly  suggest  that  the 
measured  spectral  indices  are  genuinely  mapping  the  in- situ  irregularity 
structures.  Such  issues  clearly  have  important  phenomenological 
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ATTN 
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ATTN 
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R.  Hendrick 
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y.  raien 
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R.  A-dars 
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/•  T  T  *  J 
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ATTN 
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G.  Price 

ATTN 

C. 

Reifinger 
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